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Performing highly precise Monte-Carlo simulations of SU(2) gluodynamics, we observe for the 
first time Abelian dominance in the confining part of the static potential in local unitary gauges 
such as the F12 gauge. We also study the flux-tube profile between the quark and antiquark in these 
local unitary gauges and find a clear signal of the dual Meissner effect. The Abelian electric field 
is found to be squeezed into a flux tube by the monopole supercurrent. This feature is the same 
as that observed in the non-local maximally Abelian gauge. These results suggest that the Abelian 
confinement scenario is gauge independent. Observing the important role of space- like monopoles 
in the Polyakov gauge also indicates that the monopoles defined on the lattice do not necessarily 
correspond to those proposed by 't Hooft in the context of Abelian projection. 
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Quark confinement phenomenon remains an impor- 
tant unsolved problem in quantum chromodynamics 
(QCD) ,l|. One of the most intriguing conjecture for 
its mechanism is that the QCD vacuum behaves as a 
dual superconductor due to magnetic monopole conden- 
sation 0, i.e., the color flux between a quark and 
an antiquark is squeezed into a stringlike tube as the 
Abrikosov vortex [J, [j| through the dual Meissner ef- 
fect, which yields a linear-confining potential. Although 
it is not straightforward to identify the corresponding 
monopoles in QCD in contrast to SUSY QCD Q or the 
Georgi-Glashow model @,Q with scalar fields, it is possi- 
ble to reduce SU(3) QCD into an Abelian [U(l)] 2 theory 
with magnetic monopoles by a partial gauge fixing, also 
referred to as the Abelian projection [9j, and to accom- 
modate the above dual superconductor scenario [lol ]. 

However, there arc infinite ways of the partial gauge- 
fixing. Numerically, an Abelian projection with non-local 
gauges such as the maximally Abelian (MA) gauge [ll|, 
LLJ, LL3| has been found to sup por t the Abelian confine- 
ment scenario beautifully [lj, IHj dj E3- On the other 
hand, the Abelian confinement mechanism has not been 
observed clearly so far for years in other general gauges 
in particular, in local unitary gauges [l||, [ljj [2(| . This is 
very unsatisfactory, since the quark confinement mecha- 
nism should not depend on a special gauge choice plj . 

It is the purpose of this letter to show for the first time 
that the Abelian confinement mechanism is observed nu- 
merically also in local unitary gauges with the method of 
highly precise numerical simulations. For numerical sim- 
plicity we adopt SU(2) group instead of SU(3), but the es- 
sential feature of non- Abelian gauge theory should be the 
same. As local unitary gauges, we adopt simplest candi- 
dates, namely the F12, the F123 and the spatial Polyakov 
loop (SPL) gauges as well as the original Polyakov (PL) 
gauge. Applying the multi-level noise reduction method 
invented by Liischer and Weisz (22j, we investigate the 
Abelian static potential with high accuracy and find a 



clear signal of Abelian dominance in its confining part 
in all local unitary gauges considered. Note that F12 
(F123) gauge and PL (SPL) gauge are most typical but 
are of completely different types. Since we obtain the 
same results in these different unitary gauges, we expect 
that the same can be seen also in other local unitary 
gauges. In addition, we study the flux-tube profile be- 
tween the quark and antiquark with the vacuum ensem- 
ble composed of as many as 4000 thermalized configu- 
ratio ns g enerated by means of the improved Iwasaki ac- 
tion [23[ and observe the squeezing of Abelian electric 
field into a flux tube due to the magnetic monopole cur- 
rent. This is the dual Meissner effect and the feature is 
quite similar to that already observed in the MA gauge. 
The authors expect that the results obtained here are 
very interesting to general readers, since they strongly 
suggest that the Abelian dual Meissner effect caused by 
Abelian monopoles is gauge independent and a correct 
confinement mechanism. 

We generate thermalized gluon configurations using 




FIG. 1: Abelian static potentials in the F123, the F12 and 
the SPL gauges in comparison with the non- Abelian potential. 
The solid lines denote the best fitting curves to the function 
Vm(R). 
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FIG. 2: The force in the F123, the F12 and the SPL gauges 
in comparison with the non-Abelian one. The filled (open) 
symbols are forces which are defined by using the backward 
(mid-point) difference. 

the standard SU(2) Wilson gauge action at the coupling 
constant (3 = 2.5 on the lattice TV 4 = 24 4 , where the lat- 
tice spacing is a((3) = 0.0836(8) fm. The scale is set from 
the string tension ^fa = 440 MeV. Periodic boundary 
conditions are imposed in all directions. Then we per- 
form a partial gauge fixing diagonalizing a plaquette op- 
erator X(s) = U\2(s) [F12 gauge], or an operator X(s) = 
Ui{s)U 2 {s + l)U 3 (s + i + 2)Ul(s + 3 + 2)ul{s + 3)U f 3 (s) 
[F123 gauge], or a space-like Polyakov loop X(s) = 

' p sv(2)J2l=iYlfj 1 u k( s + **) I SPL g au ge]> or a usual 
time-like Polyakov loop X(s) = U *( s + I PL 

gauge]. After gauge fixing, we decompose SU(2) link 
variables as 



U^s) = U°Js)+iaU»(s) = C^s) 





ie M (s) 



where a is the Pauli matrix, and extract Abelian link 
variables M (s) = arctan(l^( s )/t/°(s)) []]|. Since 
the first three gauges contain only space-like link vari- 
ables, we may apply the multi-level algorithm [22j to 
evaluate the Abelian static potential from the correla- 
tor of the Abelian Polyakov loop defined by Pa(s) = 



exp 



N-l 



i(s + ii) 



For the multi-level method, we 

choose the temporal extent of a sublattice to be 4a. 

Note that the expectation value of such an U(l) invari- 
ant Abelian quantity in some gauges is expressed theoret- 
ically by a sum of complicated gauge-invariant quantities 
composed of operators in various representations [24j and 
hence the numerical results are not predictable. 

In Fig. [TJ we show the Abelian static potential as well 
as the non-Abelian potential as a function of the q-q dis- 
tance R, where a tree-level perturbative improvement of 
the distance is applied to avoid an enhancement of lat- 
tice artifacts especially at short distances [25|, Hfjj]. We 
find that the results in local unitary gauges are remark- 
ably clean. Note that although the number of configu- 
rations (iV con f) for measurements is low, the multi-level 
algorithm with enough internal updates gives us clear ex- 
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FIG. 3: The profile of the Abelian electric field from the con- 
nected and disconnected correlators in the MA gauge. 

pectation values of the Abelian Polyakov-loop correlators 
up to as precisely as 10 -16 . The slope at large distances 
looks identical. We may fit the potentials to a usual 
functional form Vfn(R) = crR — c/R + [i and extract the 
string tension a. The best fitting parameters are sum- 
marized in Table |H In the FI2 and the FI23 gauges, the 
string tensions arc in almost complete agreement with 
the non-Abelian one. In the SPL gauge (iV con f = 10 with 
-^Viupd = 300000), however, the Coulombic coefficient be- 
comes so large (c w 0.8) that we cannot determine the 
string tension definitely on this small lattice. On the 
other hand, in all gauges the force, which is defined by 
differentiating the potential with respect to R, shows a 
good agreement at large distances (see Fig. [2]). In the PL 
gauge, the agreement of Abelian and non-Abelian string 
tensions is trivial, since non-Abelian Polyakov lines are 
equal to Abelian ones in this gauge. 

Next let us study the dual Meissner effect in local 
unitary gauges. In this calculation, we employ the im- 
proved Iwasaki gauge action with the coupling constant 
[i = 1.20, which corresponds to the lattice spacing a(/3) = 
0.0792(2) fm [27| . The lattice size is 32 4 with periodic 
boundary conditions. We have taken 4000 thermalized 
configurations for measurements. To improve a signal-to- 
noise ratio, the APE smearing technique is applied to the 
Wilson loop [28[. In addition, although we could measure 
usual (disconnected) correlators after the gauge-fixing, we 
evaluate a connected correlator defined by (Oa(t))w — 
(Tr[LW(R,T)tfa 3 3 A (r)])/{Tr[W(R,T)]) for various 
operators O a composed of Abelian link variables in order 
to get a better signal, where L is a product of non-Abelian 
link variables called the Schwinger line, connecting the 
Wilson loop W with the Abelian operator [2^, [3(| . r is 

TABLE I: The string tension a, the Coulombic coefficient c, 
and the constant /i obtained by the best fit. FR means the 
fitting range before tree- level improvement, x 2 is defined by 
the diagonal components of the covariance matrix and A^df 
is the number of degrees of freedom. Mupd is the number of 
internal updates used in the multi-level method. AT con f is 8. 
The errors are estimated by the jackknife method. 

aa 2 c fia FR(R/a) xV^df iViupd 



NA 
F123 
F12 



.0348(7) .243(6) 0.607(4) 3 - 10 0.35 15000 
.0350(2) .239(1) 1.187(1) 2 - 10 0.10 80000 
.0345(6) .244(4) 1.192(3) 2 - 10 1.08 80000 
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FIG. 4: The profile of the Abelian electric field in the F123. 
The solid line denotes the fitting curve for (E z (r)}w to the 
function f(r). 

TABLE II: The penetration length A and the coherence length 
£ both in unit of fm (ft = A/£), which are evaluated from 
W(5a,5a). MA(d) means the disconnected correlator in the 
MA gauge. In each case \ 2 /N^i = 0.6 — 1.1. 

MA(d) F123 F12 SPL PL 



A 



0.129(2) 0.133(3) 0.132(3) 0.134(7) 0.132(4) 
0.154(6) — — 0.162(38) 0.142(28) 



1.18(6) 



1.17(34) 1.31(29) 



the minimal distance between W and Oa ■ As Abelian 
operators, for instance, we employ the Abelian electric 
field Ei(s) = 9a{s) and the monopole current fcj(s) = 
(l/4n)e iup<T d u 9 pa (s + i) where O^s) = M (s) + 6 v (s + 
p) - M (s + v)- 9 v {s) = V(s) + 27rn^(s) (|0 M „| < tt) . 

We first examine the consistency between the result 
from the usual disconnected correlators [l6| and the con- 
nected one in the MA gauge for (E z (r))w The pro- 
file of the electric field is plotted in Fig. [31 We find 
that while the absolute magnitude of the electric field 
depends on the type of the correlator, their exponential 
decay rates look the same. Assuming a functional form 
f(r) = cicxp(— r/A) + cq, we may estimate the pene- 
tration length A, which characterizes the strength of the 
dual Meissner effect. The result is A = 0.133(4) fm for 
the disconnected correlator and A = 0.131(10) fm for the 
connected one in the case of W(5a,7a). cq is consistent 
with zero in each case. This indicates that the result of 
the connected correlator can be consistent with that of 
disconnected one. 

In Fig. 01 we show the Abelian electric field profile 
in the F123 gauge for W(5a, 5a). We find that only 
(E z (r))\y exhibits an exponential decay as a function 
of r and the penetration length is then found to be 
A = 0.133(3) fm. As seen from Table [HI A in other 
unitary gauges are almost the same, which are also con- 
sistent with that in the MA gauge [3l|. Note that we 
also investigate the profile of the magnetic field with the 
operator Bi(s) = (l/2)eijk0jk(s) and find no correlation 
with the Wilson loop. 

To identify what squeezes the Abelian electric field, 
let us study the Abelian (dual) Ampere law V x E = 
dtB + 2nk . In Fig. O we show the profile of each term 
in the F123 and the PL gauges, where only the non- 



vanishing azimuthal components are plotted. We find 
that the curl of electric field V x E is reproduced by the 
monopole current 27rfc, while the magnetic displacement 
current dtB gives only small contribution. The magni- 
tude of the profile depends on the gauge, but the qualita- 
tive feature is quite similar in both gauges. Note that this 
behavior is consistent with that in the MA gauge (2?| . 

It is important to note that the space-like monopole 
current is responsible for squeezing the electric field in 
the PL gauge. This suggests that the monopoles defined 
on the lattice do not necessarily correspond to 't Hooft's 
Abelian monopoles 0|, where the latter is due to the 
degeneracy points of eigenvalues of some adjoint opera- 
tors and they are always time- like in the PL gauge [32| . 
Rather, if the monopoles we observe on the lattice do not 
exactly correspond to the 't Hooft monopoles, the role of 
monopoles for the Abelian confinement mechanism can 
be gauge independent, and indeed, the above results seem 
to support such an expectation. 

Finally, let us briefly discuss the type of the dual su- 
perconductor vacuum. For this purpose it is impor- 
tant to evaluate the coherence length £ as well as the 
penetration length A. The ratio of these two length 
scales, the GL parameter k = A/£, classifies the vac- 
uum type. V2k = 1 corresponds to the border between 
the t ype I and the type II vacua. As demonstrated in 
Ref. [3l[ we may extract £ by fitting the profile of the 
squared monopole density (k 2 {r))w to a functional form 
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FIG. 5: Curl of the Abelian electric field, the magnetic dis- 
placement current and the monopole current in the F123 (up- 
per) and PL (lower) gauges. 
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FIG. 6: The profile of the squared monopole density (k^(r))w 
in the PL gauge for W(5a, 5a). The solid line denotes the 
fitting curve to the function g(r). 

g(r) = c[ exp(— -\/2V/£) + c o- For the operator the con- 
nected correlator is reduced to the disconnected one and 
here we evaluate the latter. In Fig. we show the pro- 
file in the PL gauge for W(5a,5a). We obtain a similar 
behavior in the SPL gauge. However, we cannot identify 
the profile in the F12 and the F123 gauges within statis- 
tics, which is probably due to contamination from many 
ultraviolet monopoles. In Table HH we summarize the 



coherence length in the SPL and the PL gauges as well 
as that in the MA gauge for comparison. The value of 
the GL parameter looks consistent with each other. How- 
ever, we note that the value may not be definite since the 
Wilson loop size adopted here is still small and further 
systematic studies are required. What we can conclude 
here is that the vacuum type of SU(2) gluodynamics is 
not far from the border. 

In conclusion, we have observed a clear signal of 
Abelian dominance in the confining part of the static 
potential in local unitary gauges for the first time. The 
structure of the flux-tube profile in these gauges strongly 
support a gauge-independent role of monopoles in the 
Abelian confinement scenario. The study of Abelian con- 
finement mechanism without performing any gauge fixing 
is in progress and the result will be published elsewhere. 
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